Testing immediate dosage compensation in Drosophila miranda via irradiation with heavy-ion beams.

Many organisms with heteromorphic sex chromosomes possess a mechanism of dosage compensation (DC) in which X-linked genes are upregulated in males to mitigate the dosage imbalance between sexes and between chromosomes. However, how quickly the DC is established during evolution remains unknown. In this study, by irradiating Drosophila miranda male flies, which carry young sex chromosomes (the so-called neo-sex chromosomes), with heavy-ion beams, we induced deletions in the neo-Y chromosome to mimic the condition of Y-chromosome degeneration, in which functional neo-Y-linked genes are nonfunctionalized; furthermore, we tested whether their neo-X-linked gametologs were immediately upregulated. Because the males that received 2-Gy iron-ion beam irradiation exhibited lower fertility, we sequenced the genomes and transcriptomes of six F1 males derived from these males. Our pipeline identified 82 neo-Y-linked genes in which deletions were predicted in the F1 males. Only three of them showed a one-to-one gametologous relationship with the neo-X-linked genes. The candidate deletions in these three genes occurred in UTRs and did not seriously affect their expression levels. These observations indirectly suggest that DC was unlikely to have operated on the neo-X-linked genes immediately after the pseudogenization of their neo-Y-linked gametologs in D. miranda. Therefore, the dosage imbalance caused by deletions in the neo-Y-linked genes without paralogs may not have effectively been compensated, and individuals with such deletions could have exhibited lethality. Future studies on sex chromosomes at different ages will further reveal the relationship between the age of sex chromosomes and the stringency of DC.

Shared evolutionary trajectories of three independent neo-sex chromosomes in Drosophila.

Dosage compensation (DC) on the X Chromosome counteracts the deleterious effects of gene loss on the Y Chromosome. However, DC is not efficient if the X Chromosome also degenerates. This indeed occurs in Drosophila miranda, in which both the neo-Y and the neo-X are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that emerged independently in D. albomicans and D. americana and reanalyzed neo-sex chromosome evolution in D. miranda Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on the neo-X is also accelerated in D. albomicans and D. americana although to a lesser extent than in D. miranda In males, neo-X-linked genes whose neo-Y-linked homologs are pseudogenized tended to be up-regulated more than those whose neo-Y-linked homologs remain functional. Moreover, genes under strong functional constraint and genes highly expressed in the testis tended to remain functional on the neo-X and neo-Y, respectively. Focusing on the D. miranda and D. albomicans neo-sex chromosomes that emerged independently from the same autosome, we further found that the same genes tend to become pseudogenized in parallel on the neo-Y. These genes include Idgf6 and JhI-26, which may be unnecessary or even harmful in males. Our results indicate that neo-sex chromosomes in Drosophila share a common evolutionary trajectory after their emergence, which may prevent sex chromosomes from being an evolutionary dead end.

Elucidating the functional evolution of heat sensors among Xenopus species adapted to different thermal niches by ancestral sequence reconstruction.

Ambient temperature fluctuations are detected via the thermosensory system which allows animals to seek preferable thermal conditions or escape from harmful temperatures. Evolutionary changes in thermal perception have thus potentially played crucial roles in niche selection. The genus Xenopus (clawed frog) is suitable for investigating the relationship between thermal perception and niche selection due to their diverse latitudinal and altitudinal distributions. Here we performed comparative analyses of the neuronal heat sensors TRPV1 and TRPA1 among closely related Xenopus species (X. borealis, X. muelleri, X. laevis, and X. tropicalis) to elucidate their functional evolution and to assess whether their functional differences correlate with thermal niche selection among the species. Comparison of TRPV1 among four extant Xenopus species and reconstruction of the ancestral TRPV1 revealed that TRPV1 responses to repeated heat stimulation were specifically altered in the lineage leading to X. tropicalis which inhabits warmer niches. Moreover, the thermal sensitivity of TRPA1 was lower in X. tropicalis than the other species, although the thermal sensitivity of TRPV1 and TRPA1 was not always lower in species that inhabit warmer niches than the species inhabit cooler niches. However, a clear correlation was found in species differences in TRPA1 activity. Heat-evoked activity of TRPA1 in X. borealis and X. laevis, which are adapted to cooler niches, was significantly higher than in X. tropicalis and X. muelleri which are adapted to warmer niches. These findings suggest that the functional properties of heat sensors changed during Xenopus evolution, potentially altering the preferred temperature ranges among species.

Transcriptomic changes with increasing algal symbiont reveal the detailed process underlying establishment of coral-algal symbiosis.

To clarify the establishment process of coral-algal symbiotic relationships, coral transcriptome changes during increasing algal symbiont densities were examined in juvenile corals following inoculation with the algae Symbiodinium goreaui (clade C) and S. trenchii (clade D), and comparison of their transcriptomes with aposymbiotic corals by RNA-sequencing. Since Symbiodinium clades C and D showed very different rates of density increase, comparisons were made of early onsets of both symbionts, revealing that the host behaved differently for each. RNA-sequencing showed that the number of differentially-expressed genes in corals colonized by clade D increased ca. two-fold from 10 to 20 days, whereas corals with clade C showed unremarkable changes consistent with a slow rate of density increase. The data revealed dynamic metabolic changes in symbiotic corals. In addition, the endocytosis pathway was also upregulated, while lysosomal digestive enzymes and the immune system tended to be downregulated as the density of clade D algae increased. The present dataset provides an enormous number of candidate symbiosis-related molecules that exhibit the detailed process by which coral-algal endosymbiosis is established.

Gene-by-Gene or Localized Dosage Compensation on the Neo-X Chromosome in Drosophila miranda.

Many organisms have a global mechanism for dosage compensation (DC) operating along the entire male X chromosome, which equalizes gene expression on the male X with that on the two Xs in females and/or on autosomes. At the initial stage of sex chromosome evolution, however, gene-by-gene (or localized) DC may also be necessary because the degeneration of Y-linked genes occurs independently at different times. We therefore tested whether the up-regulation of X-linked genes depends on the status of their Y-linked homologs, using the young sex chromosomes, neo-X and neo-Y, in Drosophila miranda. In support of the presence of gene-by-gene DC, the extent of up-regulation in males was indeed higher for neo-X-linked genes with pseudogenized neo-Y-linked homologs than for neo-X-linked genes with functional neo-Y-linked homologs. Further molecular evolutionary analysis also supports the idea that many individual neo-X-linked genes first acquired the potential for up-regulation, which then enabled the pseudogenization of neo-Y-linked homologs, without serious deleterious effects on male fitness.

Molecular phylogeny of ten intertidal hermit crabs of the genus Pagurus inferred from multiple mitochondrial genes, with special emphasis on the evolutionary relationship of Pagurus lanuginosus and Pagurus maculosus.

The hermit crab genus Pagurus exhibits high species diversity and a wide geological distribution. Despite the high species diversity of hermit crabs in the western Pacific coast of Japan, molecular phylogenetic studies of these species have yet to be conducted. To investigate their molecular phylogeny and genetic diversity, we obtained nearly complete mitochondrial genome sequences for ten Pagurus species found along the Pacific coast of Japan by next-generation sequencing, which were compared to other congeners deposited in the GenBank database. The genomes ranged from 13,458 to 16,401 base pairs in length, possessing 13 protein-coding genes, 2 rRNAs, and 22 tRNA genes. Based on the reconstructed phylogeny, we found that (1) Japanese Pagurus species separated into three groups, nested within the Northern Pacific species. (2) Pagurus lanuginosus and Pagurus maculosus, showed high morphological similarities, implying close kinship. Indeed, these two species were genetically closest to each other, compared to the remaining species studied. (3) An unspecified specimen sampled from the deep sea, which morphologically resembled Pagurus, might be a member of the Pagurus genus, but is genetically distant from the other Japanese Pagurus species. The novel data reported here may provide new perspectives for systematic studies of hermit crabs; these results provide important information that will facilitate population-level research and identifying intraspecific variation of these non-model, but ecologically important, decapod species.

Transcriptional Activation of Glycogen Catabolism and the Oxidative Pentose Phosphate Pathway by NrrA Facilitates Cell Survival Under Nitrogen Starvation in the Cyanobacterium Synechococcus sp. Strain PCC 7002.

Cyanobacteria respond to nitrogen deprivation by changing cellular metabolism. Glycogen is accumulated within cells to assimilate excess carbon and energy during nitrogen starvation, and inhibition of glycogen synthesis results in impaired nitrogen response and decreased ability to survive. In spite of glycogen accumulation, genes related to glycogen catabolism are up-regulated by nitrogen deprivation. In this study, we found that glycogen catabolism was also involved in acclimation to nitrogen deprivation in the cyanobacterium Synechococcus sp. PCC 7002. The glgP2 gene, encoding glycogen phosphorylase, was induced by nitrogen deprivation, and its expression was regulated by the nitrogen-regulated response regulator A (NrrA), which is a highly conserved transcriptional regulator in cyanobacteria. Activation of glycogen phosphorylase under nitrogen-deprived conditions was abolished by disruption of the nrrA gene, and survival of the nrrA mutant declined. In addition, a glgP2 mutant was highly susceptible to nitrogen starvation. NrrA also regulated expression of the tal-zwf-opcA operon, encoding enzymes of the oxidative pentose phosphate (OPP) pathway, and inactivation of glucose-6-phosphate dehydrogenase, the first enzyme of the OPP pathway, decreased the ability to survive under nitrogen starvation. It was concluded that NrrA facilitates cell survival by activating glycogen degradation and the OPP pathway under nitrogen-deprived conditions.

Functional divergence of duplicate genes several million years after gene duplication in Arabidopsis.

Lineage-specific duplicated genes likely contribute to the phenotypic divergence in closely related species. However, neither the frequency of duplication events nor the degree of selection pressures immediately after gene duplication is clear in the speciation process. Here, using Illumina DNA-sequencing reads from Arabidopsis halleri, which has multiple closely related species with high-quality genome assemblies (A. thaliana and A. lyrata), we succeeded in generating orthologous gene groups in Brassicaceae. The duplication frequency of retained genes in the Arabidopsis lineage was ∼10 times higher than the duplication frequency inferred by comparative genomics of Arabidopsis, poplar, rice and moss (Physcomitrella patens). The difference of duplication frequencies can be explained by a rapid decay of anciently duplicated genes. To examine the degree of selection pressure on genes duplicated in either the A. halleri-lyrata or the A. halleri lineage, we examined positive and purifying selection in the A. halleri-lyrata and A. halleri lineages throughout the ratios of nonsynonymous to synonymous substitution rates (KA/KS). Duplicate genes tended to have a higher proportion of positive selection compared with non-duplicated genes. Interestingly, we found that functional divergence of duplicated genes was accelerated several million years after gene duplication compared with immediately after gene duplication.

Genome of the pitcher plant Cephalotus reveals genetic changes associated with carnivory.

Carnivorous plants exploit animals as a nutritional source and have inspired long-standing questions about the origin and evolution of carnivory-related traits. To investigate the molecular bases of carnivory, we sequenced the genome of the heterophyllous pitcher plant Cephalotus follicularis, in which we succeeded in regulating the developmental switch between carnivorous and non-carnivorous leaves. Transcriptome comparison of the two leaf types and gene repertoire analysis identified genetic changes associated with prey attraction, capture, digestion and nutrient absorption. Analysis of digestive fluid proteins from C. follicularis and three other carnivorous plants with independent carnivorous origins revealed repeated co-options of stress-responsive protein lineages coupled with convergent amino acid substitutions to acquire digestive physiology. These results imply constraints on the available routes to evolve plant carnivory.

Evolutionary history of the extant amphioxus lineage with shallow-branching diversification.

Amphioxus or lancelets have been regarded as a key animal in understanding the origin of vertebrates. However, the evolutionary history within this lineage remains unexplored. As the amphioxus lineage has likely been separated from other chordates for a very long time and displays a marked left-right asymmetry, its evolutionary history is potentially helpful in better understanding chordate and vertebrate origins. We studied the phylogenetic relationships within the extant amphioxus lineage based on mitochondrial genomes incorporating new Asymmetron and Epigonichthys populations, and based on previously reported nuclear transcriptomes. The resulting tree patterns are consistent, showing the Asymmetron clade diverging first, followed by the Epigonichthys and Branchiostoma clades splitting. Divergence time estimates based on nuclear transcriptomes with vertebrate calibrations support a shallow diversification of the extant amphioxus lineage in the Tertiary. These estimates fit well with the closure of seaways between oceans by continental drift, ocean currents, and present geographical distributions, and suggest a long cryptic history from the origin of amphioxus to its most recent diversification. Deduced character polarities based on phylogenetic analyses suggest that the common ancestor of the extant amphioxus existed in a tiny epibenthic state with larva-like appearance of extant amphioxus, likely with ciliate epidermis.

Comparative analysis of the brain transcriptome in a hyper-aggressive fruit fly, Drosophila prolongata.

Aggressive behavior is observed in many animals, but its intensity differs between species. In a model animal of genetics, Drosophila melanogaster, genetic basis of aggressive behavior has been studied intensively, including transcriptome analyses to identify genes whose expression level was associated with intra-species variation in aggressiveness. However, whether these genes are also involved in the evolution of aggressiveness among different species has not been examined. In this study, we performed de novo transcriptome analysis in the brain of Drosophila prolongata to identify genes associated with the evolution of aggressiveness. Males of D. prolongata were hyper-aggressive compared with closely related species. Comparison of the brain transcriptomes identified 21 differentially expressed genes in males of D. prolongata. They did not overlap with the list of aggression-related genes identified in D. melanogaster, suggesting that genes involved in the evolution of aggressiveness were independent of those associated with the intra-species variation in aggressiveness in Drosophila. Although females of D. prolongata were not aggressive as the males, expression levels of the 21 genes identified in this study were more similar between sexes than between species.

Accelerated pseudogenization on the neo-X chromosome in Drosophila miranda.

Y chromosomes often degenerate via the accumulation of pseudogenes and transposable elements. By contrast, little is known about X-chromosome degeneration. Here we compare the pseudogenization process between genes on the neo-sex chromosomes in Drosophila miranda and their autosomal orthologues in closely related species. The pseudogenization rate on the neo-X is much lower than the rate on the neo-Y, but appears to be higher than the rate on the orthologous autosome in D. pseudoobscura. Genes under less functional constraint and/or genes with male-biased expression tend to become pseudogenes on the neo-X, indicating the accumulation of slightly deleterious mutations and the feminization of the neo-X. We also find a weak trend that the genes with female-benefit/male-detriment effects identified in D. melanogaster are pseudogenized on the neo-X, implying the masculinization of the neo-X. These observations suggest that both X and Y chromosomes can degenerate due to a complex suite of evolutionary forces.

Two-step evolution of endosymbiosis between hydra and algae.

In the Hydra vulgaris group, only 2 of the 25 strains in the collection of the National Institute of Genetics in Japan currently show endosymbiosis with green algae. However, whether the other non-symbiotic strains also have the potential to harbor algae remains unknown. The endosymbiotic potential of non-symbiotic strains that can harbor algae may have been acquired before or during divergence of the strains. With the aim of understanding the evolutionary process of endosymbiosis in the H. vulgaris group, we examined the endosymbiotic potential of non-symbiotic strains of the H. vulgaris group by artificially introducing endosymbiotic algae. We found that 12 of the 23 non-symbiotic strains were able to harbor the algae until reaching the grand-offspring through the asexual reproduction by budding. Moreover, a phylogenetic analysis of mitochondrial genome sequences showed that all the strains with endosymbiotic potential grouped into a single cluster (cluster γ). This cluster contained two strains (J7 and J10) that currently harbor algae; however, these strains were not the closest relatives. These results suggest that evolution of endosymbiosis occurred in two steps; first, endosymbiotic potential was gained once in the ancestor of the cluster γ lineage; second, strains J7 and J10 obtained algae independently after the divergence of the strains. By demonstrating the evolution of the endosymbiotic potential in non-symbiotic H. vulgaris group strains, we have clearly distinguished two evolutionary steps. The step-by-step evolutionary process provides significant insight into the evolution of endosymbiosis in cnidarians.

Different Endosymbiotic Interactions in Two Hydra Species Reflect the Evolutionary History of Endosymbiosis.

Endosymbiosis is an important evolutionary event for organisms, and there is widespread interest in understanding the evolution of endosymbiosis establishment. Hydra is one of the most suitable organisms for studying the evolution of endosymbiosis. Within the genus Hydra, H. viridissima and H. vulgaris show endosymbiosis with green algae. Previous studies suggested that the endosymbiosis in H. vulgaris took place much more recently than that in H. viridissima, noting that the establishment of the interaction between H. vulgaris and its algae is not as stable as in H. viridissima. To investigate the on-going process of endosymbiosis, we first compared growth and tolerance to starvation in symbiotic and aposymbiotic polyps of both species. The results revealed that symbiotic H. viridissima had a higher growth rate and greater tolerance to starvation than aposymbiotic polyps. By contrast, growth of symbiotic H. vulgaris was identical to that of aposymbiotic polyps, and symbiotic H. vulgaris was less tolerant to starvation. Moreover, our gene expression analysis showed a pattern of differential gene expression in H. viridissima similar to that in other endosymbiotically established organisms, and contrary to that observed in H. vulgaris We also showed that H. viridissima could cope with oxidative stress that caused damage, such as cell death, in H. vulgaris These observations support the idea that oxidative stress related genes play an important role in the on-going process of endosymbiosis evolution. The different evolutionary stages of endosymbiosis studied here provide a deeper insight into the evolutionary processes occurring toward a stable endosymbiosis.

Evolutionary Transitions of MicroRNA-Target Pairs.

How newly generated microRNA (miRNA) genes are integrated into gene regulatory networks during evolution is fundamental in understanding the molecular and evolutionary bases of robustness and plasticity in gene regulation. A recent model proposed that after the birth of a miRNA, the miRNA is generally integrated into the network by decreasing the number of target genes during evolution. However, this decreasing model remains to be carefully examined by considering in vivo conditions. In this study, we therefore compared the number of target genes among miRNAs with different ages, combining experiments with bioinformatics predictions. First, we focused on three Drosophila miRNAs with different ages. As a result, we found that an older miRNA has a greater number of target genes than a younger miRNA, suggesting the increasing number of targets for each miRNA during evolution (increasing model). To further confirm our results, we also predicted all target genes for all miRNAs in D. melanogaster, considering co-expression of miRNAs and mRNAs in vivo The results obtained also do not support the decreasing model but are reasonably consistent with the increasing model of miRNA-target pairs. Furthermore, our large-scale analyses of currently available experimental data of miRNA-target pairs also showed a weak but the same trend in humans. These results indicate that the current decreasing model of miRNA-target pairs should be reconsidered and the increasing model may be more appropriate to explain the evolutionary transitions of miRNA-target pairs in many organisms.

Experimental Approach Reveals the Role of alx1 in the Evolution of the Echinoderm Larval Skeleton.

Over the course of evolution, the acquisition of novel structures has ultimately led to wide variation in morphology among extant multicellular organisms. Thus, the origins of genetic systems for new morphological structures are a subject of great interest in evolutionary biology. The larval skeleton is a novel structure acquired in some echinoderm lineages via the activation of the adult skeletogenic machinery. Previously, VEGF signaling was suggested to have played an important role in the acquisition of the larval skeleton. In the present study, we compared expression patterns of Alx genes among echinoderm classes to further explore the factors involved in the acquisition of a larval skeleton. We found that the alx1 gene, originally described as crucial for sea urchin skeletogenesis, may have also played an essential role in the evolution of the larval skeleton. Unlike those echinoderms that have a larval skeleton, we found that alx1 of starfish was barely expressed in early larvae that have no skeleton. When alx1 overexpression was induced via injection of alx1 mRNA into starfish eggs, the expression patterns of certain genes, including those possibly involved in skeletogenesis, were altered. This suggested that a portion of the skeletogenic program was induced solely by alx1. However, we observed no obvious external phenotype or skeleton. We concluded that alx1 was necessary but not sufficient for the acquisition of the larval skeleton, which, in fact, requires several genetic events. Based on these results, we discuss how the larval expression of alx1 contributed to the acquisition of the larval skeleton in the putative ancestral lineage of echinoderms.

Microsatellite marker development by multiplex ion torrent PGM sequencing: a case study of the endangered Odorrana narina complex of frogs.

The endangered Ryukyu tip-nosed frog Odorrana narina and its related species, Odorrana amamiensis, Odorrana supranarina, and Odorrana utsunomiyaorum, belong to the family Ranidae and are endemically distributed in Okinawa (O. narina), Amami and Tokunoshima (O. amamiensis), and Ishigaki and Iriomote (O. supranarina and O. utsunomiyaorum) Islands. Because of varying distribution patterns, this species complex is an intrinsic model for speciation and adaptation. For effective conservation and molecular ecological studies, further genetic information is needed. For rapid, cost-effective development of several microsatellite markers for these and 2 other species, we used next-generation sequencing technology of Ion Torrent PGM™. Distribution patterns of repeat motifs of microsatellite loci in these modern frog species (Neobatrachia) were similarly skewed. We isolated and characterized 20 new microsatellite loci of O. narina and validated cross-amplification in the three-related species. Seventeen, 16, and 13 loci were cross-amplified in O. amamiensis, O. supranarina, and O. utsunomiyaorum, respectively, reflecting close genetic relationships between them. Mean number of alleles and expected heterozygosity of newly isolated loci varied depending on the size of each inhabited island. Our findings suggested the suitability of Ion Torrent PGM™ for microsatellite marker development. The new markers developed for the O. narina complex will be applicable in conservation genetics and molecular ecological studies.

Development and characterization of 27 new microsatellite markers for the Indian bullfrog Hoplobatrachus tigerinus and its congeneric species.

The Indian bullfrog Hoplobatrachus tigerinus and its four congeneric species are common frog species and distributed throughout South Asia. Due to recent human activity, they are facing a changing environment and reduction in natural population size. For effective conservation and molecular ecological studies, we therefore isolated and characterized microsatellite loci for these frogs. We obtained genomic data using an Ion Torrent PGM sequencer and designed 54 primer sets for candidate loci. By screening for polymorphic loci in individuals of H. tigerinus and its congeneric species, we isolated 27 loci as highly polymorphic microsatellite loci. Eight of these loci were commonly applicable for all species except H. chinensis. Within two populations of H. tigerinus, the total number of alleles per locus and expected heterozygosity ranged from 2 to 18 and 0.271 to 0.938, respectively. No significant linkage disequilibrium was observed across all loci, and five showed a significant deviation from Hardy-Weinberg equilibrium in some populations after Bonferroni correction. Consequently, our findings suggest that these novel markers will be applicable for conservation genetic studies across varying scales from inter-population to inter-individual.

Development and characterization of 14 microsatellite markers for Buergeria japonica (Amphibia, Anura, Rhacophoridae).

Buergeria japonica is a common frog species distributed throughout almost all islands in Ryukyu Archipelago. Because of their exceptionally wide distribution and higher physiological tolerance comparing to the other anurans, their demographic history and formation of distribution are intrinsic topics in the herpetological fauna of Ryukyu. Microsatellite marker is ideal genetic marker for such studies at inter- and intra-population level. We therefore developed microsatellite markers of B. japonica utilizing Ion PGM™ sequencing. As a result of the screening, we developed a total of 14 polymorphic markers. To test availabilities of these markers, we genotyped four island populations. The total number of alleles and expected hetelozygosities per locus ranged from 4 to 21 and 0.00 to 0.864, respectively. The phylogenetic relationship among the four populations based on the genetic distances of these markers was congruent with general divergence pattern of amphibians and reptiles in Ryukyu area. These markers developed in this study are considered to be useful for future studies about phylogeography and demography of this species.

Tissue- and stage-dependent dosage compensation on the neo-X chromosome in Drosophila pseudoobscura.

Sex chromosome dosage compensation (DC) is widely accepted in various organisms. This concept is mostly supported by comparisons of gene expression between chromosomes and between sexes. However, genes on the X chromosome and autosomes are mostly not homologous, and the average gene expression level on these chromosomes may not be the same even under DC, which complicates comparisons between chromosomes. Many genes with sex-biased expression also make comparisons between sexes difficult. To overcome these issues, we investigated DC by comparing the expression of neo-X-linked genes in Drosophila pseudoobscura with those of their autosomal orthologs in other Drosophila species. The ratio of the former to the latter in males would be 1 under DC, whereas it becomes 0.5 without DC. We found that the ratio was ∼0.85 for adult whole bodies, indicating that the DC is incomplete on the neo-X chromosome in adults as a whole. The ratio (∼0.90) was also significantly less than 1 for adult bodies without gonads, whereas it was ∼1.0 for adult heads. These results indicate that DC varies among tissues. Our sliding-window analysis of the ratio also revealed that the upregulation of neo-X-linked genes in males occurred chromosome wide in all tissues analyzed, indicating global upregulation mechanisms. However, we found that gene functions also affected the levels of DC. Furthermore, most of the genes recently moved to the X were already under DC at the larval stage but not at the adult stage. These results suggest that DC in Drosophila species operates in a tissue/stage-dependent manner.